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We study the HBT interferometry of ultra-relativistic nuclear collisions using a freezeout model 
in which free pions emerge in the course of the last binary collisions in the hadron gas. We show 
that the HBT correlators of both identical and non-identical pions change with respect to the case of 
independent pion production. Practical consequences for the design of the event generator with the 
built in Bose- Einstein correlations are discussed. We argue that the scheme of inclusive measurement 
of the HBT correlation function does not require the symmetrization of the multi-pion transition 
00 amplitudes (wave-functions). 
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ryy , I. INTRODUCTION 

-xt- : 

Pion interferometry is expected to provide important information about the space-time picture of ultrarelativistic 
J> ' nuclear collisions. It has already proved to be a sensitive tool for the detection of the collective motion of the matter 
created in the course of ultrarelativistic nuclear collisions [|l]-[| . The primary goal of this paper is to show that under 
CN| ' favorable circumstances, interferometry is capable of detecting the difference between various mechanisms of freezeout 
' which is the last transient phase before the regime of free streaming of hadrons. The effect we rely upon is entirely due 
to real interactions at the freezeout stage, and its magnitude may serve as a measure of these interactions. We predict 
that the normalized two-particle correlator of non-identical pions (and even for pairs like irp, pn, etc.) must differ 
fS from the reference unit value at any difference Afc of the pions momenta. The normalized correlator of identical pions 
O ' does not approach unity at large Afc and slightly exceeds value of 2 at Afc = 0. (In the idealized model which is used 
in this paper to demonstrate these effects, they are small.) The second goal is to discuss the effect of the multiparticle 
final states on the one- and two-pion inclusive spectra. We demonstrate that the contribution of these states depends 
on the microscopic dynamical mechanism of freezeout. We explicitly show that the effect of the multi-particle final 
states is limited by the actual range of the interactions at the freezeout stage. No particles, which are causally or 



(and) dynamically disconnected from the two-pion inclusive probe, can contribute this effect. Our general conclusion 
is that the problem of the analysis of HBT data cannot even be posed without reference to an explicit dynamical 
model. Finally, our model may serve as a prototype of a realistic dynamical mechanism of freezeout, provided that 
the preceding stage of evolution is a hot gas of hadrons. Our results can be used for the simulation of truly quantum 
two-pion distributions with the input from the event generators based on semi-classical dynamics. 

The paper is organized as follows: In Sec. [nj the effect is explained at the introductory level. We pose the problem 
of interferometry (in its most rigorous form) as a problem of a quantum transition specified by the observables and 



the initial data in Sec. III. Sec. IV deals with the formulation of the pion interferometry problem and the results 
of calculations for the hydrodynamic theory of multiple production. A more realistic model of freezeout, which 
incorporates some elements of hadron kinetics, is considered in Sec. The analytic answer for the model-based 



calculations is obtained by the end of this section. In Sec. VI, we discuss the effect of multi-pion final states on the 



inclusive two-pion spectrum and speculate about physical phenomena that may lead to the increase of the effect. A 



practical issue of the event generator, with the built-in Bose-Einstein correlations, is discussed in Sec. VII 
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II. PHYSICAL MOTIVATION 



The simplest and most popular introductory explanation of the principles of HBT interferometry is as follows: The 
amplitude to emit ("prepare") the pion at the point xpj and to detect it with the momentum k\ is ajv 
For a system of two pions prepared at points xn and xm 1 and detected with the momenta k\ and k 2 , the transition 
amplitude (wave function) is a superposition of two indistinguishable amplitudes, 

A NM {ki,k 2 ) = fl^^Je-^^ttMtfeJe-* 1 " + aArfeJe^^aMftje-^'" , (2.1) 

and for a system of distributed pion sources, the two-particle inclusive spectrum is 

' LX '" } =2 V {| aw (fc 1 )a M (fc2)| 2 +Re [(aAr(fc 1 )a M (fc 2 )a^(fc2)a; f (fc 1 ) e ^ fc ^ fc2 )( a: "- a; -)]} . (2.2) 



dkidk-2 



N.M 



This scheme carries an implicit assumption that the pions are truly independently created in the state of free 
propagation. Unfortunately, in the literature, this property is mostly taken for granted, though in the real world, it 
takes place only in very special occasions. For example, independence of pion production in nuclear collisions has 
been used to describe Bose-Einstein correlations at Bevalac energies when the excited nuclear matter is dominated 
by nucleons In that case, the bremsstrahlung of pions, that accompanies scattering of nucleons, was suggested to 
provide independent pion sources. 

At RHIC energies (100 GeV/nucleon), we anticipate that the nuclear matter before the freezeout is totally dominated 
by pions and only very few nucleons are expected in the central rapidity region. There is a long-standing conjecture 
that the last form of the matter, before the freezeout, is a hot expanding hadronic gas. If this is indeed the case, then 
the main type of interaction at the freezeout stage is the binary collisions of pions. Therefore, the free propagation of 
any pion starts after it has experienced the "last" collision in the hadronic gas, and there is no really independent free 
pion production. Indeed, in this case, at least four particles are involved in the interference process, because at least 
two particles appear in the final state in each binary collision. If two pions (say, ■k + tt + ) are detected, then there are 
at least two more undetected particles emerging from two different collisions, and these particles may be identical as 
well. Thus, we encounter an additional (hidden) interference which affects the measured inclusive two-particle cross 
section.^] Moreover, even if the detected particles are different (e.g., ir + ir~), their partners still may be identical (e.g., 
the process ir + ir° — > ir + n° takes place at the coordinate xjv, and 7r~7r° — > ir~ir a takes place at the coordinate xm)- 
Therefore, if the last interaction is the binary collision in the pion gas, the naive scheme (2.1) has to be modified. In 
this modified case, there are four interfering amplitudes, 

A NM (fci , fc 2 ; qi , g 2 ) = a N (h , q x )a M {k 2 , 9 2 )e~ l(fcl +qi)xN e -«(«*+»)*" 
+a N {k u q 2 )a M {k2,qi)e- i{kl+q2 > N e- i ^ + ^ XM 
+a N {k 2 ,q 1 )a M {k 1 ,q2)e- l ^ + ^ XN e-^ + ^ XM 
+a N (k 2 ,q 2 )a M (k uqi )e-^ + ^ XN e-^ + ^ x ™ . (2.3) 



1 We essentially base this on the Dirac definition of the wave function as a transition amplitude, ipk(x) = (k\x) = tp*(k), which 
allows for two complementary interpretations: (i) the particle is prepared with momentum k and detected at the point x, or, 
equivalently, (ii) the particle is prepared at the point x and detected with the momentum k. The same interpretation remains 
valid for the entire hierarchy of the multi-particle wave functions, {ki, kt\xi, £2), etc. 

2 A similar mechanism is known to provide corrections to the collision term in kinetic equations as well as the first virial 
corrections to the equation of state of quantum gases pi . 
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dk\dk 2 



The expression for the two-particle spectrum thus becomes more complicated, 

X] \ANM(ki,k 2 ;qi 1 q 2 )\ 2 = 4 ^ ^ { \a N (ki, qi)a M (k 2 , g 2 )| 2 

JV,M gi,92 AT,M <ji.?2 

+Rc [ a^(fci,«i)o^(fc2, gaJo^Cfea, g x )a^ {k u g^e^-M^*™) 
+a N (k 1 ,q 1 )aM(k2,q2)a%(k u q 2 )al J (k 2 ,q 1 )e- i ^-^ XN - x ^ 
+ aw (A :i ,g 1 )aM(fc2,g2)aw(fc2,g2)aM(fci,'7i)e- l(fel+9l - fc2 -' ?2)( ^- a;M) ] } , 



(2.4) 



and contains three interference terms (see Fig. [l]). The first interference term corresponds to the familiar "opened" 
interference in the subsystem of the detected pions. It is present only if the detected pions are identical. The two 
other terms describe the "hidden" interference in the subsystem of undetected pions which, however, are unavoidable 
partners of the detected pions in the process of their creation in the final states. These terms are present even if the 
detected pions are different. In fact, we deal with the dynamically generated situation when the system of two pions 
is described by a density matrix and not by the wave function. When the dynamics of the system is driven by binary 
collisions, the two pions just cannot be found in a pure stateQ 
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FIG. 1. Diagrams for four terms in Eq. (2.4). Bold solid and dashed lines correspond to the detected and non-detected pions 
in the final state, respectively. The second graph is absent when the detected pions are non-identical. The vertical dotted lines 
cut the diagrams through the on-mass-shell lines corresponding to the final-state pions. The external legs correspond to the 
initial-state pions. 

A full scenario of ultrarelativistic heavy-ion collisions is still absent and an understanding of the quantum kinet- 
ics from microscopic quantum field theory, rather than from semi-classical approach, is only beginning to emerge. 
Currently, the most self-consistent scenario is based on an assumption that equilibration happens very early. Then 
it becomes possible to appeal to the equation of state of nuclear matter at various stages of the scenario, starting 
from the quark-gluon plasma (QGP) and ending up with the hadronic gas (as is done, e.g., in Ref. 0.) If this is 
not the case (no hadronic gas occurs), the pions may be created in the course of hadronization and immediately 
freely propagate. In this situation, it is more likely that the pions are created independently. Thus, the difference 



between Eqs. (2.1), (2.2) and (2.3), (2.4) becomes of practical importance. Pion interferometry provides a tool which 
is capable of distinguishing between these two scenaria. In Sec. jv|, we argue that the effect of hidden interference 
may become even stronger if the system approaches the freezeout stage with the "soft" equation of state when the 
correlation length is long (as is, e.g., in the vicinity of a phase transition). 

Furthermore, it is easy to understand that the correlations due to the hidden interference exist even, e.g., between 
pions and protons, protons and neutrons, etc. All these effects, which may carry significant information about the 



At high pion multiplicity, we can neglect an exceptional case when two detected pions come from the same binary collision. 
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freeze-out dynamics, require (and deserve) special study, which is beyond the scope of this paper. We insist that it 
is highly desirable (though not easy) either to measure all these correlation effects or to establish their absence at a 
high level of confidence. Both results are important since they would allow one to constrain the full self-consistent 
scenario of heavy-ion collisions which is absent now. These constraints may also affect theoretical predictions of the 
photon and dilepton yields in heavy-ion collisions. 

III. THEORETICAL BACKGROUND: DEFINITION OF OBSERVABLES 

A precise definition of observables is extremely important, because the HBT interferometry does not allow one to 
pose the mathematically unambiguous inverse problem. In general, we have to formulate a model and compare the 
solution with the data, relying on common sense and physical intuition. The full set of assumptions that accompany 
the formulation of the model is never articulated in full. Particularly, the question about the nature of states in which 
the particles are created is never discussed. However, this is the key issue. Interferometry is a consequence of the 
interference, and, in quantum mechanics, the latter cannot be even addressed without direct reference to the quantum 
states. Indeed, interference takes place every time when, with a given initial state, there are at least two alternative 
histories of evolution to a given final state. The wave function is nothing but a transition amplitude between the 
two states. HBT interferometry studies the two-particle wave functions (or, more precisely, the two-particle density 
matrix). Therefore, the most natural way to avoid any ambiguity is to pose the whole problem as a problem of a 
quantum transition.^ 

Let | in) be one of the possible initial states of the system emitting a pion field which has three isospin components 

tt(x) ee {TTi(x)} = {ir + (x),ir (x),ir_(x)} , (3.1) 

where n^(x) — tt^_(x), and no(x) = tt})(x). At the moment of the measurement (tf — > +00) each component of the 
pion field can be decomposed into the system of analyzer eigenfunctions f k {x) , 

tt 4 (x) = J dk\Ai(k)f h (x) +A\{k)f* k {x% h(x) = (2Tr)-^ 2 (2k )- 1 / 2 e- tk - x . (3.2) 

This expansion holds only after freezeout. In general, tt(x) is a multiplet of Heisenberg operators driven by the 
evolution operator S. The annihilation operators Aj(fe) of the pion field are given by 

Mk) = f d?xfUx) idl 7Ti(x) , A\(k) = f d 3 xirl(x) id* f k (x) . (3.3) 

Jx°=tf Jx°=tf 

The operator Ai(k) describes the effect of a detector (analyzer) far from the point of emission, so, by definition, the 
pion is detected on mass-shell, k° = (k 2 + to 2 ) 1 / 2 . The inclusive amplitudes to find one pion with momentum k and 
two pions with the momenta k\ and ki in the final state are 

(X\Ai(k)S\m) and (X|A i (fe 1 ) J 4 J -(fe 2 )5'|in) , (3.4) 

respectively. Here, the states \X) form a complete set of all possible secondaries. Summing the squared moduli of 
these amplitudes over all (undetected) states \X), and averaging over the initial ensemble, we find the one-particle 
inclusive spectrum 



4 We closely follow (especially, in Sec. [V]) Ref. g, where this approach was first used for photon interferometry of the QGP. 
There, the photon emission was due to the processes qq — > 717 and qg — > jq, and the identical gluons or quarks of the final state 
would interfere as well, if they had not become the constituents of the thermalized system. 
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= Trp itl y2^A\(k)\X)(X\Mk)S = Trp in ^Al(k)Mk)S , (3.5) 
dk ^ 

where the density operator /5;„ describes the emitting system. Of these two equations, the second one describes the 
algorithm of the measurement. In other words, we deal with the operator A^(fc) = A\(k)Ai(k), which gives the 
number of pions of the z'th kind detected by an analyzer tuned to momentum k. The first equation indicates that 
all effects of multiparticle production are accounted for in this basic definition of the inclusive one-particle spectrum. 
The multiparticle states are encoded in a sum over the complete set of the unobserved states, 

^2\X){X\ = l, (3.6) 
x 

where they enter with an a priori assumption that they are identically weighted, and all together form a unit operator. 
In the same way, one may obtain the inclusive two-pion spectrum 

r(2) 



= Trp in Y / S^Al(k 1 )A}(k 2 )\X)(X\A J (k 2 )A l (k 1 )S = TTp^ . (3.7) 



Once again, this equation defines the number of pairs, 

NaikiM) = A\(k 1 )A](k 2 )A j (k 2 )A i (k 1 ) = JVi(fci)[jy,(fc 2 ) - <M( fe i - fc 2)] , (3.8) 

as the observable, and incorporates all multiparticle effects by its derivation. This is a function of actual dynamics 
which is driven by the evolution operator S to select the states which physically contribute to the process of the 
measurement, and to assign them the dynamically generated weights. 



Equations (3J>) and (3/7) are universal in a sense that the standard observables of interferometry are expressed in 
terms of their Hcisenberg operators. In order to make them useful, we have to specify both the evolution operator S 
and the initial data embodied in the density matrix p ln . In other words, the theory that claims to have any predictive 
power must incorporate physical information about the freezeout dynamics. Two models which reflect our current 
vision of the heavy-ion scenario, are explored in the next two sections. 



IV. NAYVE FREEZEOUT. 



The simplest and the most naive model of freezeout has been explained in detail in Ref. Q. In order to have a 
reference point for the discussion of the more involved situation detailed in the next section, we review this model 
below with new updated emphases. In fact, this model is not dynamical. No microscopic mechanism of free pion 
production is specified. We just declare the pion field to become free after some time (or, in a more intelligent 
way, starting from some space-like surface S c ). This is what is done technically, but certain assumptions must be 
kept in mind. First, up to the freezeout surface S c , the system is assumed to be a continuous medium which obeys 
relativistic hydrodynamic equations. The hypersurface E c corresponds to some critical temperature T c at which all 
interactions that maintained the local thermal equilibrium in the expanding medium are switched off. Consequently, 
we assume the pion distribution over the momenta (not in phase space!) at this moment to be (almost) thermal. It 
is incorporated into the initial data for the future free propagation. 

Second, the correlation length along S c is finite and much less than the total size of the system which is a dynamical 
consequence of the interactions before the freezeout. In order to incorporate this property into the naive picture of 
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the instantaneous freezeout, we must consider the whole system as a collection of boxes (fluid elements) filled by free 
particles which are opened when their world lines reach E c . 

Third, even though we may wish to disregard the quantum nature of the pions in the fluid expansion phase, we 
have to account for it when we pose the problem of interferometry. In other words, the pions have to be produced in 
certain states. This (perhaps, most important) goal is also achieved if we mimic freezeout by the model of the opening 
boxes. 



Practically, we act as follows: The field Tti(x) in Eqs.(3.3) has to be evolved in time starting from the initial data 
on the hypersurface E c , i.e., 

n^x) = [ dE„(y) G rct (x - y) & m(y) . (4.1) 



The space of states in which the density matrix acts is also defined on this surface. Substituting Eqs. (4.1) to (p.q), 
we find the pion Fock operators, expressed in terms of the initial fields, 

S^A(k) S = Jd 3 x f* k (x) JdZ^y) G mt (x - y) % m(y) . (4.2) 

This equation may be simplified using the explicit form of the free pion propagator, 

&Mk)S = 9(x° - y°) [ d^(y)f* k (y) id% wi(y) . (4.3) 



The answer is simple because the whole problem of evolution is reduced to just the free propagation of the pion field. 
Now, we have to incorporate the idea of the absence of long-range order on E c . This is done in three steps. First, we 
replace the continuous integral over E c by the sum of integrals over the "cells" labeled by a discrete index N, 



dZ(y) - V / d 3 y , (4.4) 

Sc M JV* 



where Vf is the three-dimensional volume of the iV-th cell on the hypersurface E c . (Hereafter, the ^-labeled quantities 
are related to the local reference frame with the time-axis normal to the hypersurface E c .) The pion field within the 
7V-th cell also acquires an additional label N . Next, we perform the second quantization of the pion field within each 
cell independently, 

TJv(y) = ^2[a N (p)(t) p {y) + a\{p)(j)* p {y)\ , 
p 

4> p (y) - (2 V£por 1/2 e- <p - a . Pl=P 2 + ™ 2 , (4.5) 

where the components of vector p take discrete values defined by the boundary conditions on the walls of each box, 
and the ^-labeled quantities are related to the local rest-frame of a fluid element. Finally, we impose the commutation 
relations on the pion field, 

[a N {p 1 ),a\ I {p 2 )] = 8 NM S Pl p 2 , (4.6) 

which is equivalent to the dynamical independence of the fields belonging to different cells. Each pion is created and 
propagates independently of all others. Thus, the model is completely defined. It is mathematically very simple and 
reflects the main features of the physical process. There are several scales in this model, the size I of the elementary 
cell, the freezeout temperature, T c , and the geometric parameters L of the flow. (The correlation length is of the 
order I > 1/T C and, in general, we must require that I <C L.) The interplay of these parameters should be explicitly 
accounted for in the course of calculations. The one-particle spectrum of pions is 
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At this point, we must treat one of the integrals as a delta- function which sets the momenta p& and k^ equal, while 
the second integral becomes just the volume Vff of the fluid cell. This procedure requires that |fe|, \p\ » l/£ ~ T c , 
i.e., the measured momenta should be sufficiently high. In the same way, we obtain the expression for the two-particle 



spectrum. Since by virtue of Eq. (4.6) we have 

(a N (pj alj (p 2 )a M > (p' 2 (Pi ) ) = 
= {Snn'Smm'S PiP ' 1 5 P2P > 2 + d N M'SMN'Sp lP ' 2 S P2P 'J(a N (p 1 )a' N (p 1 )}(aM(P2) a M(P2)}, 

the two-pion inclusive spectrum becomes, 

dN^f dNW rfiV (1) , ,>),,,, / \ t / \ \ 

" - 2^ 2^ (MPi)<r(Pi)>( a M(p 2 )<r(P2)) 



dk\dk 2 dk\ dk 2 

NM pi ,P2 



cos(fci — k 2 )(xN — xm) 



# 



d 3 ye- l{kl - pi)y / d 3 ye +l(fcl " P2)y / d 3 ye' l{k2 ' P2)v / d 3 ye +l{k2 - pi)v . (4.8) 



Thus, we have reproduced Eq. ( |2.2| ) which followed from an intuitive conjecture (2.1) about the interference of two 
indistinguishable amplitudes. These amplitudes, by their design, reflect all relevant properties of the interaction which 
prepare the two-pion system and a device which detects this system. This single device consists of two detectors (e.g., 
two tracks) tuned to the pions with momenta k\ and k 2 . Since, by the definition of the inclusive measurement, nothing 
else is measured, there are two and only two interfering amplitudes. Regardless of how large the total number, A/", 
of pions produced in a particular event is, only the two-pion transition amplitude has to be symmetrized since the 
remaining M — 2 pions are not measured. 

In order to avoid any misunderstanding, we remind the reader that the quantum-mechanical measurement of some 
observable, by definition, includes the procedure of averaging over an ensemble. A single element of this ensemble 
carries no quantum- mechanical information. Following Ref. ||, we can argue that the inclusive measurement explores 
all quantum mechanical fluctuations which can dynamically develop before the moment of measurement and are 
consistent with the detector response. From this point of view, nothing but a causal chain of real interactions can 
affect the detector response, and only dynamical histories can interfere. There is a significant difference between the 
cases when one-, two-, or three-particle inclusive distributions are measured. Each of these measurements is unique 
in a sense that they are all mutually exclusive even if they are obtained from the same ensemble of multiparticle 
events. For example, if the three-pion distribution is measured (i.e., the average over an ensemble is accomplished), 
then the integration over the momentum of the third pion does not result in the two-pion inclusive distribution p2[ . 
This example reflects a qualitative difference between classical and quantum distributions. The former are defined 
immediately in terms of probabilities while the later are defined via transition amplitudes with the probabilities playing 
the secondary role. 

The same conclusions follow from the field-theory formalism we employ. The interference emerges as a strict 
consequence of the commutation relations (JAM) which incorporate a distinctive property of the hydrodynamic model 



to localize the freezeout point for each pion independently. We start from the quantum operator (3.8) of the measured 



observable and trace the two-pion signal back to its origin. Eventually, we arrive at the square of the symmetrized 
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two-pion amplitude, once again, regardless of the total number A/"of pions in a particular event. The states with many 



pions are completely accounted for in Eqs. ( p. 4 ) and (3/7) and, consequently, in (4/7) and {<L8). However, they can 



contribute to the one- and two-pion observables only via real interactions which are absent in our oversimplified model 
of a naive freezeout. A more realistic example with the interaction is considered in the next section. 



In Eq.(4.8), we seemingly encounter a problem. If, e.g., the first of the integrals is considered as the delta- function 
which sets momenta pf and kf equal, then the value of the second one is not obvious. It occurs that the second 
integral over the volume becomes 

d 3 ye +l{kl - k2)v = V*. (4.9) 

To prove this, we need more physical information. The interferometry would have been impossible if we could have 
traced each pion back to the coordinate of its emission or, in other words, if we could have built an "optical image" 
of the source. According to the Rayleigh criterion, the latter is possible only if |fci — k?\ L > 1 . Therefore, we may 
observe interference only provided the last inequality does not hold. Thus we have |fci — k 2 \ < 1/L <C which 



proves ( |4.9| ) 



In the continuous limit, by introducing an auxiliary "emission function" J(k\,k 



2) 



J(h,k 2 ) = I dE^x)^^-n(ki-u(x))e- i ^- h »> , (4.10) 



the expressions for the one- and two-particle inclusive spectra can be rewritten in a form which is convenient for 
numerical computations: 

, m 



dk 



J{k,k) , (4.11) 



and 



fc?fcg dN \ = J(fci,fci) J(fc 2 ,fc 2 ) + Re 



dk\dk 



2 



J (hi, fe) J(&2, k\ 



(4.12) 



respectively. 

In Ref. [||, these equations were used to study interferometry for several types of one-dimensional flow. Here, we 
are interested only in the case of the boost-invariant geometry which must be very close to the reality of RHIC. Indeed, 
any scenario initiated with a strong Lorentz contraction (up to 0.1 fm!) of the nuclei cannot possess a scale associated 
with the initial state. Hence, both at classical and quantum levels, the system must evolve with the preserved boost- 
invariance. In terms of the hydrodynamic theory of multiple production Jic| , the absence of scale in initial data of the 
relativistic hydrodynamic equations immediately leads to the Bjorken self-similar solution p"l| ] as the only possible 
solution. Recent analysis of quantum fluctuations at the earliest stage of heavy-ion collisions ^ indicates the same. 
The overlap of Lorentz-contracted nuclei converts them into a system of modes of expanding plasma with the global 
boost-invariant geometry. This is a single quantum transition and it is not compatible with the picture of gradually 
developing parton cascade |1^| . 

Theoretical analysis of two-pion correlations, in the picture of "naive freezeout," from the phase of expanding hot 
pion gas was done in Refs. fjlH- One of the predictions of these studies, the so-called m^scaling of the pion and kaon 
correlators, was confirmed by the data of the NA35 and NA44 collaborations obtained from Au-Au collisions 
at SPS energies (^10 GeV/nucleon). These data carry two important messages. First, even when the nuclei are 
Lorentz-contracted only up to the size of 1 fm, this is almost enough to bring about the boost invariant regime of 
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collective flow. There is no doubts that this picture will be even more pronounced at RHIC. The second, less trivial 
consequence of the observed mt-scaling is that the freezeout is sharp; if the creation of the final-state pions were 
extended in (local) time, the m t -scaling would vanish p4[ . Thus, addressing the pion production at RHIC, we have 
every reason to rely on two facts: (i) the hadronic matter before the freezeout forms a collective system, and (ii) this 
system is in the state of the self-similar boost-invariant expansion. 

For the immediate goals of this study, an advantage of the intensive longitudinal flow is that it provides a window 
in the phase-space of the two pions where the effect of the hidden interference (explained in Sec. and computed 
in more detail in Sec. [v|) is most visible. This window corresponds to the measurement of the transverse size of the 
longitudinally-expanding pipe using pairs of pions with the same rapidity. In this way, we employ a distinctive feature 
of the hydrodynamic-type sources to localize the emission spectrum. The required localization in rapidity is achieved 
by choosing pions with large transverse momenta (p t > 3T ~ Zm^). This has a very simple physical explanation. If 
the emitting system is kinetically equilibrated (or even sufficiently chaotized) then the mean energy per particle is 
limited from above. Thus, in the rest-frame of a fluid element, the particles with large transverse momenta can have, 
on average, only small longitudinal momenta. In other words, such particles are effectively frozen into the collective 
hydrodynamic motion in the laboratory frame. 5 Measuring there longitudinal rapidity after freezeout, we are most 
likely to measure the longitudinal rapidity of the fluid at the emission site. 

The parameters of the model with the Bjorken geometry are the critical temperature, T c (~ m^), and the space- 
like freeze-out hypersurface, defined by t 2 — z 2 = t 2 = const. The coordinates and the four-velocity of the fluid are 
parameterized as x M = (t cosh 7/, r, r sinhry), and it^x) = (cosh 77, 0, sinh?7), respectively. The rapidity, 77, of a fluid cell 
is restricted to ±Y in the center-of-mass frame. We assume an axially-symmetric distribution of hot matter in a pipe 
with area S± = txR\. The particles are described by their momenta, k% = (k®,ki,k!-) = (mi cosh 9i , ki , sinh 9i ) , 
where ki is the two-dimensional vector of the transverse momentum, 9i is the particle rapidity in z-dircction, and 
m| = m^ = m 2 + kf is the transverse mass. Let us introduce 2a = 61 — 62, 26 = 61+62- The one-particle distribution 
is expected to be close to a thermal distribution of pions at temperature T = T C . Since the m± values of interest are 
larger than T c , we may take this distribution in a Boltzmann form, and use the saddle-point method to estimate the 



integrals (4.12). For the one-particle distribution these two steps yield 



dN 



dB\dk\ J-y 
The general expression for J(fci, k 2 ) is 



S±mi / drj cosh(6q — rj)e 



mi cosh(0i-7j)/T c 



rS±mi 



2nT c 
m 1 



-mi/T c 



(4.13) 



J(fci, k 2 ) 



i(k 1 -k 2 )r 1 



T I drj 

-v 



(mi + ma) cosh(# — 77) cosh a — (mi — 7712) sinh(77 — 6) sinh a 



x exp 



1 

fi. 



(mi + iF(ni\ — 1112)) cosh(?y — 6) cosh a — (mi + iF(m\ + 7712)) sinh(7; — 6) sinhi 



(4.14) 



where F — tT c . Once we consider only those pions with large transverse momenta, the integral over the rapidity 77 
can (and should) be computed in the saddle-point approximation. This yields 



R(k lt k 2 ) 



dN 2 



dy 1 dp 1 dy 2 dp 2 
x exp 



dTVi dNx 



dyidpi dy 2 dp 2 



h(z) cos 



1 - 1 f(\ff IP ^ 9{z)g{l/z) 
H{z) 1. H(l/z) 1 



5 This is the physical origin of the so-called mt-scaling and the qualitative basis for the saddle-point calculations below. 







' , , . H(z) , ,1, . H(l/z) 
h(z) sin + h(-) sin v / ; 



+ G(z) + G(-) 



(4.15) 



where fx = (miTOa) and z = (mi^) 1 ' 2 , and we have introduced the functions 



h(z) = 



tan H (z) 



~F 2 (z- i) 2 + 4F 2 sinh 2 a 



-4f 2 (z 2 



(z 2 + cosh 2a) 2 + F 2 (z 



cosh 2a)' 
1 , 



1/4 



1/2 



2F(cosh2a 



-F 2 (z- ±) 2 +4.F 2 sinh 2 a 



tan G(z) 



F(z 2 - 1/z 2 ) 
z 2 + cosh 2a 



(4.16) 



When the transverse momenta of two pions are equal, then mi = 7712, and Eq. (LIE) reproduces the result of Ref. [Q. 
The function / depends on how we parameterize the transverse distribution of the matter in the longitudinally 
expanding pipe. If f(r) = 8{R 2 L - r 2 ) then f(\Ak x \R±) = [2 Ji(|AA;| j R l )/|AA;| j Rl] 2 , where J x is a Bessel function. 
This HBT correlator possesses a natural property of any two-point function in the boost-invariant geometry; it depends 
on the difference of rapidities and not on the difference of the longitudinal momenta. The mt-scaling of the "visible 
longitudinal size" is just a synonym for this property, which supports a picture of the sharp freeze-out of matter in 
the state dominated by the boost-invariant-like longitudinal expansion. 
Two remarks are in order: 

1. The shape of the correlator given by Eqs. (4.15) and ( 4.16| ) is manifestly non-Gaussian. Thus, it is plausible to 
get rid of the intermediate Gaussian fit and to use these equations to fit the data as the first step. 

2. A main deficiency of the correlator (4.15) is that the possible transverse expansion has not been taken into 



account. This is not easy to do since the initial data in the transverse plane are, as yet, poorly understood and the 
amount of matter involved in the collective transverse motion at the freezeout stage may vary depending on the initial 
conditions, equation of state, etc. We consider this issue still open. 

V. DYNAMICAL FREEZEOUT. 



Now, let us formulate a more realistic model for the freezeout, relying on the following phenomenological input. 
Let the system before the collision be a pion-dominated hadronic gas. Therefore, since pions are the lightest hadrons, 
the regime of continuous medium is supported mainly due to 7T7r-scattering. In the framework of chiral perturbation 
theory there are two main channels contributing to this process. The first channel is s-wave scattering which we shall 
model using the sigma-model. The second channel is the p-wave 7T7r-scattering via the p-meson. We shall use a model 
interaction Lagrangian of the form 



(5.1) 



where cr(x) is the field of the scalar a-mesons, and p M = (p~ , p° , p + ) is the field of the p-mesons. The reader can 
easily recognize here that part of the Lagrangian of the sigma-model which is necessary to reproduce the tree-level 
(or even skeleton) amplitudes of the s- and p-wave 7T7r-scattering. The evolution operator is usually defined as 



S = Texp{i J C int (x)d i x} 



(5.2) 



Applying the standard commutation formulae ul 
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A^S-SA^k) = J d^x-^f^x) , SU^fc)-^)^ = | d V fc (^^^y , (5.3) 



to the Eqs. (3J3) and (3/7) we obtain 

dfc J \<57r [a] (yi) J7ij o] (y 2 

for the one-particle spectrum of the pions of kind a (ir + , 7r°, 7t~), and 

diVPA r / <5 2 5 t c5 2 5 \ 

V = / d 4 y 1 d i y 2 d 4 y 3 d 4: y4fk 1 (yi)fk 2 (y3)( j j—rj 1 — v ~r~: r~ r~ , )/£, (ife)/fc 2 (y4) > ( 5 - 5 ) 

dfe l dfe 2 7 \<JT[o](yi)*7T[6](y3) ^(ya)^ (y 4 )/ 

for the spectrum for pairs of pions of kinds a and b. According to our agreement, the angular brackets denote an 
average weighted by the density matrix p m . [In order to avoid confusion, we place the isospin indices in square 
brackets.] 

Since the goal of this paper is to demonstrate a physical effect, we shall limit ourselves, in what follows, with the 
it interaction term. In this case, we are able to do most calculations analytically, which is necessary in order to 
understand the interplay of the many parameters. More involved calculations will require a realistic description of 
the 7T7r-scattering and Monte-Carlo computation of multiple integrals. 

Calculation of the one-pion spectrum is relatively simple. A direct computation of the functional derivatives in 
Eq.(|J) leads to 

[a] _ Hiir I ,4 j4„. e 



/Jk y -f / rfWy 2fc V(Tt(^t . n^(x)^(x)n^(x) :)T(: ^\y)^\y)^\y) : 5)) , (5.6) 

where the symbols T and denote the time and anti-time orderings, respectively. Coupling the field operators 
according to the Wick theorem (to the lowest order of the perturbation expansion we have to put 5=1), the average 
in Eq. ( ^ ) becomes, 

(•••) = * 3 [4?WM^^ - (5-7) 

where any pair of arguments x and y lie within the same cell of a size defined by the correlation length in the system 
under consideration. This is a consequence of the cellular structure of the density operator p- ln which determines the 
average (• • •). The correlators GAB{x,y) are defined (using the Keldysh technique [jl6|; see also Ref. ||) as 

G%l(x,y) = -i(T c {^{x A )^\y B )) . (5.8) 

Taking for a = 1 (tt + ) and using the facts that cx 6 t: * , G^(y,x) — G^\x,y), and G^\y,x) — G^\x,y) we 
arrive at 

(•••)= i 3 G [ ^\y,x)[3G^\x,y)G^\y,x) + G [ ^(x,y)G^ ] (y,x)] . (5.9) 



Next, we have to insert this expression into Eq. (|5.6|) and pass over to the momentum representation (cell- by-cell) 
according to 

G^(p) = Gf a 3] (-p) = -2m[9(p )n^(p)+0(-p a )} , 
G 1 W ] (P) = G l £ ] (-p) = -2m{6(p )+9(-p°)n^(p)} , 
Gg? ] (p) = -2wi[6(p°)nW(p) + 6(- P )} , 
G [ ™ ] (p) = -2m[6(p°) + 0(-p> ( 0) (-p)] . (5.10) 
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These equations reflect a very important feature of the process under investigation (already accounted for in 



Eqs. ( p.4| ) and (5^5). ) Namely the initial states in the system of colliding pions are populated with the densities 
n ffl(p) °f cnar g e d pions and n$(p) of neutral pions. All the final states of free propagation are not occupied. They 
are virtually present with an a priori unit weight in the expansion ( |3.6D of the unit operator. Access to the states with 
many pions is provided dynamically by real interactions, and these states show up in higher orders of the evolution 
operator expansion. Depending on the type of commutation relations, these states always appear in a properly 
symmetrized form, and there is no need to start with any ad hoc symmetrization. By definition, the wave function is 
the transition amplitude between the prepared initial state and detected final state. Unless the transition is explicitly 
measured, the wave function does not exist at all and there is no object for symmetrization. Therefore, there is no 
statistical correlations associated with the multiple production of pions. All correlations are dynamical (see further 



discussion in Sec. |VIp 



After some manipulations, we obtain for the one-particle spectrum 



dN K+ Sivi/W f d Pl d P? d 1 X n i W 2 2 W 2 2s x , 2 2n 

~dfe~ = T^ " J 2fcQ(2^)8 «S( fc + g-Pi-P2)3(Pi-m )S(p 2 -m )6(q -m ) 

x{3n { +\p 1 )n ( J\p 2 ) + 6n { +\p 1 )n { N \p 2 ) + 2n { +\p 1 )n { °\p^ . (5.11) 



Here, as in the model of the "naive freezeout" discussed in Sec. IV, one of the integrals over the space-time volume 



of the domain (where the individual collision occurs), is treated as the delta- function (viz., momentum conservation), 
while the second integral becomes just the four-volume of the interaction domain. (More details are given in 
Appendix A.) In what follows, we shall limit ourselves to a simplified picture when all the statistical weights of the 
initial state are taken in the Boltzmann form and the chemical potential is zero. In this case, we gain an overall factor 
of 3+6+2-1-1=12 and arrive at 



dk 

N 



£ 4 4) / ^kOiS)^ ^ + 1*)*(P? - m2 WP* ~ m 2 )e-^"" (pi+P2) . (5.12) 



Now, the integration over I = p\ — pi can be explicitly carried out (see Eq. (B.l) 



«± . £«4V*"~> / ^f^w - ^ - ^ . (, 13 , 



Introducing an auxiliary function, 



I(k,N,x NM ) = J d 4 q5+(q 2 - m 2 )e-^«) e*™ Jl - j^f-^ , (5.14) 

which is computed in Appendix B, we may write explicitly 

dN$ Uvgl _ 4 4) * N(kuN) Tfu AT ns Ugl.K^m/T) 



dk 16 ^ 2fc (27r) 8 V ' ' ; 16 m/T J 2k (2w) 8 



e -ku(x)/T 



iwmM^ , ,5.5, 
16(27r) s to/2 V m ± 



where, with reference to the data discussed in Sec IV, we assume that the freezeout occurs in a very small interval 
At of the proper time r. Additionally, we further simplify the model by taking the temperature the same throughout 
the entire freezeout domain. 

Calculation of the two-pion spectrum is very cumbersome but follows exactly the same guideline, 
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dN 



(2) 



dkidk2 16 



= d^ V2 dVy A - 



-ik 1 (y 1 -y 2 ) p -ik 2 {v3-y4) 



4fc0fc^(27r) 6 

xT(:7T^(y 2 )7r^(y 2 ) ^\y 2 ): :^ 4 V tW (zmV 1 " 1 M^)) , (5.16) 
where, as before, we must put S = 1 in the lowest order of the perturbation expansion. If the observed pions are 



identical, then the result of coupling is a very long expression which is given in its full form in Appendix A, Eq. (A.l). 
The next step is to integrate over the elementary volumes Vjv- This is done according to Eq. (A.3). Finally, using 
the Boltzmann approximation for the statistical weights, we carry out an explicit integration over l\ = s\ — s 2 and 
h 



Pi — P2 (see Eq. (B.l)). The result is as follows, 

dN {2 } , dN {1 ) dN {1 } TT 2 qj v - (4) (4) a ru > fl n, ^ 

7T+7T+ _ 7T + 7T+ , Jin \ Tr (4) „ (4) -fj N ( kl U N ) „ ~PM ( &2 U M ) 

dk x dk 2 dki dk 2 + 1fi 2 ^ N 

d i qid i q 2 



N,M 



4fc0fc0(27r) 16 



4m 2 
(ki + qif 



Am 2 
(k 2 + q 2 ) : 



. e -PN(qiu N ) e -0M(q2UM) §+ ( q 2 _ m ^)S + (q 2 - m 2 ) 



xRe{144 e~ i ( kl ~ k2 ^ XN ~ XM ') + 44 e ~ i ( 9l_92 ) (xjv_:EM ) + 4 : 4 :e ~ i (ii-+ k i--<i2-k 2 )(x N -XM)^ 



(5.17) 



In this equation, we recognize the symmetrized transition amplitude for all four pions engaged in the process of 
the two pion production. Each term acquires a factor associated with the isospin algebra. Though all multi-pion 



states are included in the definitions (5.4) and (5.5) of the observables, it is neither necessary nor even possible to 
symmetrize them, if only one or two pions are measured. The next orders of the evolution operator expansion will 
bring in additional final-state particles. They will all be properly symmetrized. However, these terms will be more 
and more suppressed (discussion in Sec. |v| ). 

In a similar way, the original expression, (A. 2), is transformed to the inclusive spectrum of two different pions, 

( 2 ) jw-W jjwW „4 



dN 



dkidk 2 



dN -+ dN i- , Win \^ T/(4) 1/ (4)„-/3 N (fc 1 u JV )„-/3 M (fc 2 u M ) 

fc7~ 162 N M e 

N,M 



dk\ dk 2 



d i qid 4 q 2 
4fc0fc0(27r) 16 ' 



Am 2 



Am 2 



e -p-N(qiu„) e - P M(<l2UM) S f q 2 _ m ^)S + (q 2 - m 2 ) 

(k 2 + q2r 

xRe[16 e ~ i ( qi ~ q2 ^ XN ~ XM " 1 +26 e ~ i ( qi+kl ~ q2 ~ k2 ^ XN ~ XM 



(5.18) 



Next, we have to find the net yield of the interference terms after the momenta of the unobserved pions are integrated 
out. As a first step, we can use an auxiliary function I(k,x, Ax), ( 5.14 ) and rewrite Eqs. ( 5.17 ) and (|5.18|) , replacing 
the discrete sums by the integration over the continuous distribution of the collision points, 

r(2) 



dN y 



dk\dk 2 



1 4 

* 9%* 
16 2 



d^x\d^x\ 



e ~(k 1 u{x 1 ))/T Ce -(k 2 u{x 2 ))/T c 



4fc0fc°(27r) 16 

xRe{144/(fci, Xl , 0)J(fc 2 , x 2 , 0)[1 + e - l (M-k 2 )(x!-x 2 ) 
-44/(fe 1) Xl , (xi - x 2 ))P{k 2 ,x 2 , (xi - x 2 ))[l + e -<(*i-*»)(«i-^)]} 



(5.19) 



for identical pions, and 



dN 



(2) 



dk\dk 2 



1 A 

16 2 



d A x\d A x\ 
4fcOfcO(27r) 16 



e ~(k 1 u(x 1 ))/T Ce ~(k 2 u(x 2 ))/T c 



xRe{UAI(k 1 ,x 1 ,0)I(k 2 ,x 2 , 0) + I(ki,x u (z x - x 2 ))P (k 2 ,x 2 , (n - x 2 ))[16 + 2&e' l{kl ~ k2)(xi ~ X2) ]} 



(5.20) 
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for non-identical pions. These equations are written in the same approximation as the one-pion inclusive spectrum 



( E-15D - 

The first two terms of Eq. ( 5.1 9| ) have a well known form of the product of the intensities of two independent 
sources and an interference function [1 + cosAkAx] and correspond to the standard scheme of HBT interfcrometry. 
In these terms, the space-time integration over coordinates x\ and x 2 is factorized and reproduces the result of the 
naive model of freezeout given by Eq.(4.15). The only difference is an inessential kinematic factor and a form-factor 



I(k 1 ,x 1 ,0)I(k 2 ,x 2 ,0) = 4ir 2 



K x {m/T) 



m/T 



(5.21) 



which eventually cancels out in the normalized correlator. In the last two terms of Eq. ( [5.1S| ), we encounter an 
additional factor T{x\,x 2 ) — I{k\, X\, Ax)I*(k,2 1 x 2 , Ax) which, in general, does not allow one to factorize the space- 
time integrations. In our special geometry of the freezeout, however, some further simplifications are possible, provided 
the transverse momenta kit and km are large. Indeed, integrating with respect to xi, we end up with a modified 
version of the emission function ( 4. 10| ) , 



J(h,k 2 ;x 2 ) = / d*x 1 e- kl < Xl M T e- i < hl - k >) Si r(x 1 ,x 2 ) 



At / d 2 ri/(ri)e" ?1 ^ 1 ~*^ / Tdr)ie~^ mi ^ cos Hm- e i)- iT [ m i cosh(j7i-0i)-m 2 cosh(rji-e 2 )] F(Ar),Ax) 



-Y 



where the form-factor T can be conveniently rewritten as 



T(Arj, p) = 47r 2 m- 



K 1 (mV~U 2 ' 



(5.22) 



(5.23) 



where 



ir = — - 

y2 



4Tt{Tt + i) sinh 



2V1-V2 



and we denoted p = fi — r 2 . The saddle point of the integration over 771 in Eq. ( 5.22| ) is defined by the equation,^] 



tanh(?7i — 0) 



1 + iTr[l + (m 2 /mi)] 



tanh a 



(5.24) 



1 + iTr[l - (m 2 /mi) 

and is not affected by the form-factor ( 5.23| ). Thus, if we select pions with the same rapidity, then 2a = 6\ — 2 = 
and the saddle point occurs at the point rji = 8 = 6\ = 62- Integrating in the same way with respect to x 2 , we come 
to the conclusion that 772 = 6. The form-factor then becomes independent of rapidities 77, and its argument becomes 
a real function^] 



T{0,p) = ATr 2 m 



2 ^(mva/n+p 2 ) 



m 2 [(l/T 2 



(5.25) 



From now on, we will deal only with pions of the same rapidity, which means that we physically take aim at two 
scattering processes which occur in one narrow slice of the longitudinally expanding pipe. In this particular case, the 



very cumbersome formulae (like (4.15)) become drastically simplified. The spectrum of two identical pions reads as 



6 T he ex ponential in Eq.(3.22) is the same as in Eq.(4.14). Thus, it is exactly the saddle point (5.24), which leads to Eqs. (4.15) 
and ( |4.16D . 

7 When 771 = 772, the form-factor T is a smoothly decreasing function. When 771 7^ 772, it acquires an additional oscillatory 
pattern in the A?7-direction. 
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dN 



(2) 



dk\dk 2 
x 1 Re 



dM 1 } dN (1 } 



2 4 
T Sir 



144 



16 2 4fcO/fc°(27r) 16 



4^ 2 m 4 (Ar) 2 r 2 



2ttT 
1 

y/m 1 m 2 



3 -(mi+m 2 )/T c 



^/TOiTO 2 



[m\m 2 + F(F — i)[mi — w^) 2 ] 1 / 2 

44 /" 
+ Y^2tt / r 1 dr 1 r 2 dr 2 f(r 1 )f{r 2 ) 



2vr / rdrf(r)J (|Afc|r) 
2 - K 2 (mV(l/T 2 )+p 2 ) 



K\{m/T) 



m/T 

M\Ak\ P ) 



m 2 [{l/T 2 )+p 2 l 
2 * K 2 (m^{l/T 2 )+p 2 ) 



44 /" 

144 7 Jo rrr [(1/T 2 ) + p 2 

where p 2 = r 2 + r| — 2rir 2 cos0. In the same way, we obtain the two-particle spectrum of non-identical pions, 



(5.26) 



(2) 



dk\dk 2 



dN^dN^_ j ^gl 

dki dk 2 



144 



4tt 2 to 4 (At) 2 t 2 



2ttT 



16 2 4fc?fc0(27r) 16 

x27r f ridrir 2 dr 2 f{ri)f(r 2 ) 
16 26 



-(mi+ma)/T c 







x < 1 Re 

144 144 



^Jm\m 2 



yfmim 2 

2 ^ KUm^(l/T 2 )+p 2 ) 
m 2 [(l/T 2 )+p 2 } 

Jo(|Afc|p) 



\m\m 2 + F(F — i)(mi — TO2) 2 ] 1 / 2 UVI ltv J ^ 
In order to give a quantitative estimate of the effect (yet with all reservations intact due to the imperfection of our 
model description of the 7T7r-interaction) we shall assume that the distribution of the longitudinally expanding hot 
pipe in the transverse direction is homogeneous within a cylinder of radius R± . We also neglect the transverse flow. 
Next, we normalize the two-pion spectrum by the product of the one-pion spectra and introduce 

dN^>/dkxdk 2 



C a b(ki,k 2 )- 

(dN?> /dk^dN^ /dk 2 ) 

Finally, we consider only the case of |fci| = \k 2 \. Then, m\ — m 2 and we arrive at 

' 2Ji(|Afc|i?) \ 2 t 11 f R ri dr ir2 dr 2 f 2 ^ 2tt K 2 [(m/T)y/l + T 2 p 2 l 



\Ak\R ) 



(i 



(nR) 2 J q " K 2 (m/T) 



l + T 2 p 2 



1 



R r 1 dr 1 r 2 dr 2 f 2 * 2ir K 2 [(m/T) y/l + T 2 p 2 ] 
# —, ' " 



i + H M\Ak\ P ) 



(5.28) 

[1 + M\Ak\p)} , (5.29) 
(5.30) 



/„ (7ri?) 2 J r K 2 {m/T) l + T 2 p 2 

where we remind the reader that the longitudinal rapidities of the pions are set equal. The remaining integrations in 
Eqs. ( 5.29Q and ( 5.30D are performed numerically and the results are plotted in Figs. |^ and [| 
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FIG. 2. Normalized correlator of identical pions as a function of q — (Afc) s id c (with (Afc) out = and A6* = 0): (1) mR± = 1; 
(2) mR± = 2; (3)- mR± = 3. The right plot is a magnified portion of the left plot which allows one to see the scale of correction 
due to the hidden interference. 
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VI. POSSIBLE CONTRIBUTIONS OF THE MULTIPARTICLE STATES. 



As we have seen in the previous section, the multi-pion final states indeed contribute to the one- and two-pion 
inclusive cross sections. This contribution is due to the real interactions, and therefore, we may speculate about 
its magnitude in different situations and at different conditions of observation. Within the unrealistic 7r 4 -model we 
employed to make analytic solution possible, the effect is understandably small. Indeed, the additional pions emitted 
in the collision process lead to the form- factor ( [5.25| ), F{p), in the interference term. Let us consider the formal limit 
of T — > oo. Then the form- factor becomes 

Tf-m \ A 2 i K l{ m \fP) 

r(p, T — ► oo) = 47r to — ^- . (6.1) 

m z p z 

This expression is immediately recognized as the square of the pion vacuum correlator, 

Gx{p) = -7(0|7r(pV(0) + ^(0)7r(p)|0) , 

which represents the density of states in the pion vacuum. These are exactly the states which are excited as the final 
states of free propagation in the course of the 7T7r-collision, and the characteristic scale of correlation for these vacuum 
fluctuations can be nothing but m~ . At T — > oo, the colliding pions are (formally) very hard (their thermal wave- 
length becomes very short) and they cannot explore any scale except this one. Thus, the "hidden interference" is active 
only if two collisions occur within the range of a typical fluctuation in the "trivial" pion vacuum. Otherwise, these 
collisions are dynamically disconnected and the interference becomes strongly suppressed. At finite temperatures, the 
behavior of the form- factor at the origin is less singular, but the main scale for the distance p remains JW" 1 , which 
defines the smallness of the effect in the oversimplified model tt 4 . The parameter which regulates the magnitude of the 
excess of the normalized correlators ( 5.29| ) and ( 5.30f ) is the ratio of the effective radius of the dynamical correlations 



to the full transverse size R± of the expanding pipe. As is seen from the plot of Fig. [|, at T = the excess reaches 
0.08 when R± ~ m" 1 . 

This understanding of the nature of the scale of correlation shows that the magnitude of the effect may be quite 
different when the 7T7r-interaction is treated more adequately. For example, with properly fitted parameters, the er- 
model provides a good description of many data which depend on 7T7r-interaction. Most of low-energy interactions of 
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elementary particles are mediated via wide resonances like /o (which is sometimes identified with er), a±, p, etc. In this 
picture, the 7T7r-interaction is not as local as in the 7r 4 -model, and more scales show up resulting in the dependence 
of the 7T7r-cross section on the momentum transfer. The low-energy dynamics is rich and includes a chiral phase 
transition. In the vicinity of this phase transition the correlation radii should increase. Thus, on the one hand, we 
may speculate about various dynamical phenomena which may lead to an increase of the correlation radii. On the 
other hand, the system may develop a strong transverse hydrodynamic motion. In this case, the "measured -Rj_" 
also becomes a dynamically defined effective quantity, which is (though not as much as in the longitudinal direction) 
smaller than the natural transverse size. 

Until now, we have considered the lowest order (with respect to the interaction) contributions to the one- and 
two-pion inclusive distributions. The next orders can be accounted for by the expansion of the evolution operator in 
the Eqs. (5J:) and ( |5.16| ). Even without explicit calculations, it is clear that there will be two types of corrections. The 



virtual radiative corrections should be absorbed into the definitions of the dressed propagators and vertices. Within 
the framework of our phenomenological approach, these are of no interest or value. The real corrections are connected 
with the production of additional pions and involve integration over the positions of the additional points of emission 
(collision) and over the momenta of the additional non-observed pions. The corresponding complicated transition 
amplitudes will be properly symmetrized. However, the result of the integration over the unobserved momenta will 
again be a set of form-factors which will cut off all distances between the points of scattering which exceed the 
correlation length. Thus, the higher order effects will always be suppressed with respect to the lowest order. 

Our last remark is about the possibility of the induced emission of pions. To address this question, one has to 
realize that the so-called Bose-enhancement has two different (and mutually exclusive) aspects. If we deal with the 
measured states, then the symmetrization is solely due to the measurement. This is the essence of the HBT effect. If 
the states are not measured, then we indeed deal with the induced emission which is due to the factor \/n + 1 in the 
equation a^n) = \/n + l\n+ 1) , which requires that the process of the emission of the (n + l)-th quantum physically 
occurs in the field of all n identical quanta emitted before. This is the basic idea behind the laser. If the pions are 
produced with space-like separation, then stimulated emission is impossible. 

VII. BOSE-EINSTEIN CORRELATIONS IN EVENT GENERATORS 

Our estimates show that the effect of hidden interference is not expected to be large unless the system undergoes 
the freezeout starting from a soft mode with long-range correlations. This possibility is very appealing due to its 
possible physical richness. Careful simulations will be required in order to estimate the signal-to-noise ratio in a real 
detector. This job is usually done with the aid of the so-called event generators. A complexity stemming from an 
additional interference makes simulations more challenging. We address this issue below. 

Over the last decade, many computer codes which are supposed to mimic the particle yields in heavy-ion collisions 
have been developed. These codes rely on classical propagation of the particles in hadronic cascades. They do not 
describe the quantum properties of the measurement which lead to the interference of amplitudes. At best, the 
output of an event generator is the list of particles with one-to-one correspondence between their momenta and the 
coordinates of their emission. Such a list is equivalent to a complete optical image of the source. Therefore, there are 
no alternative histories in the free propagation. An "afterburner" has to be added to incorporate quantum effects. 

In order to recover the alternative histories and thus simulate the two-pion correlations in a model with independent 
pion production, one has to know four one-particle amplitudes, ajv(k)- They are necessary to recover the two-pion 



amplitude -4tva/ (&i, fe) (see Eq. (2.1) ). Hence, one has to store some information about the pion sources before 
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the emission which allows for such a reconstruction. In the model of binary collisions, the parameters of parents 
(their ID, momenta, and the coordinates of the reaction volume) have to be stored for each pion from the detected 
pair. Otherwise, it is impossible to reconstruct all four amplitudes which interfere in the process of the two-pion 
measurement. 

Eventually, one has to average the square of the inclusive amplitude over an ensemble of all events. Then, this 
procedure has to be repeated for every point in the two-particle momentum space. In more sophisticated models, with 
the pions coming from decays of resonances, we may have even more interfering amplitudes and each case requires 
special consideration |n| . 

The requirement that the quantum mechanics of the last interaction, which forms the two-pion spectrum, must be 
recovered, follows from first principles and cannot be removed. In some particular cases, one may take a short cut 
and weaken the requirements of the complete description. However, every precaution should be taken that the initial 
state of the interfering system is described as a quantum state. 

Finally, we have to mention that in the literature, another kind of correlation function, which has been introduced 
2p| l ) in connection with the problem of simulation of Bose- Einstein correlations in multi-pion 
events, is widely discussed. The algorithmic definition of this object (which is different from the inclusive differential 
distribution used in this paper) has several steps: (i) The events are sorted by the total number TV of pions in the 
event, (ii) The correlation function CV(fci, fc/v) is measured, i.e., the A/"-pion amplitude is fully symmetrized and its 
squared modulus is averaged over the subset of events with this given Af. (iii) All but two arguments of Cjv(fci, k_\f) 
are integrated out, which results in a set of two-pion correlators CV(fci, A^), still dependent on M. (iv) This set is 
averaged over A/", in order to obtain C(ki, fo). We could not find a quantum observable which would correspond to 
this correlator. We also could not establish a direct connection between this correlator and the parameters of the 
emitting system. 



L§ (see also Refs. 19 , 



VIII. CONCLUSIONS 



In this paper, we show that precise measurements of two-particle correlators, both for identical and non-identical 
particles, may uncover important information about the freezeout dynamics. A reliable theoretical prediction for the 
magnitude of these correlations can be made only within the framework of an elaborated scenario of the heavy-ion 
collision which is currently absent. Data may be very helpful for the theoretical design of the fully self-consistent 
scenario. 

It is equally important to measure the magnitude of the effect of hidden interference or to establish its absence 
at a high level of confidence. From this point of view, our calculations carry an important message that there are 
no theorems which require the normalized two-pion correlator to approach unity for large differences of momenta. 
Neither should it take the value of 2 at Afc = 0. The latter value is the norm (charge) of the state with two identical 
pions in a pure state only. In general, two final-state pions are described by a density matrix. For the same reason, 
the normalized correlator of non-identical pions should not equal 1. Moreover, the correlator should not equal 1 even 
for the couples like 7rp, pn, etc. Therefore, the technique of "mixing of events" should be used with extreme caution, 
if used at all. The idea of a universal Gaussian fit of the HBT correlators can hardly be fruitful as well. (The authors 
fully realize that there may be many technical difficulties in the practical implementation of this advice.) 

One more practical lesson of our analysis concerns the theory of multi-pion correlations. This issue cannot even 
be addressed without an explicit account of the dynamical mechanism responsible for the creation of multi-particle 
states. Our general conclusion is that HBT interferometry cannot be model- independent. 
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APPENDIX A. 



Here, we reproduce lengthy expressions for the two-particle inclusive spectra as they appear after performing all 
couplings. We do not use them in this form throughout the paper. However, they are the starting point, would we 
wish to design a code for the simulation of the two-pion correlations. In this code, the particle densities n^ip) will 
serve as the frequencies for the initial-state pions to appear within the range of the last interaction. At least some part 
of the integrations with respect to the momenta of the unobserved final-state pions is preferable to do analytically in 
order to simplify the remaining Monte-Carlo integrations. 

The full expression for the two-particle spectrum of identical pions as it appears after performing all couplings in 
Eq.(^16|), is 



dN 



(2) 



dN # dN$ gl v 
16 
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xe 
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where s — s\ + S2, p = Pi + P2; %n and Xm are the "central" coordinates of two cells, and xnm 
long expression emerges for the correlator of two different pions as well. 
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+ e - l ( ?1 +fc 1 - 92 -fe 2 )^M [17n ( + ) (si)n (^) (s2)n ( + ) (pi)n ( i -) (p2) 

-An i ^\si)n { ~\s 2 )n { ° I \p 1 )n { ^ (p 2 ) + An ( °\s 
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(A.2) 



In Sec. [v], we proceed with a simplified version of these equations which emerges when the distributions n N ^ (p) are 
taken in Boltzmann form with an additional assumption that the system is locally neutral. This results in a common 



weight function and coefficients which are just the numbers of terms in different groups in Eqs. (A.l) and 



Integration over the elementary volumes in Eqs. (A.l) and (|A. 2|) requires special discussion: There is an important 



difference between the volumes Vn in equations of Sec. IV and the volumes Vat here. The former are related to the 



elementary fluid cells where the distributions un(p) of particles over their momenta are established due to many 
collisions inside the iV'th cell , and we do not (and even are prohibited to) localize the coordinates of individual 
collisions within the cell explicitly. The latter correspond to the opposite case, when we locate the space-time 
coordinates of individual collisions. Therefore, the volumes Vat correspond to the actual range of the pion interaction 
potential and are much smaller than the volumes Vn of the fluid cells required by the hydrodynamic picture. Hence, 
the meaning of the distribution functions tin{p) changes as well; now they correspond to the probability for the pion 
with momentum p to penetrate the "reaction domain" near the space-time point xn- In fact, the volume Vat is 
defined by the cross section of the 7T7r-interaction itself. This picture is consistent only if the pions themselves are 
considered as wave packets and not as plane waves. The volume Vn is that volume where the incoming packets are 
identified by the interaction and where the outgoing packets are completely formed. With this picture in mind, we 
can integrate, e.g., 

d yi e~ l{kl+qi - s)vi [ dy 2 e t{kl+q2 - s]y2 = (2 7 r) 4 ,5(fc 1 + q x - s) j dy 2 e- l{qi - q2)v2 = (27r) 4 <5(fc 1 + q x - s)V N . (A.3) 

V N Jv N JV N 

In the first equation, we assume that the relation between the volume of integration and momenta allows one to verify 
the conservation of momentum in the collision. In the second equation, we took into account that the actual range 
for the momenta q x and q 2 is of the order of 2T ~ 2m n ^ while the length of the ir-ir scattering is ~ l/5m n . In all 
subsequent calculations we can rescale Vn back to the elementary volume Vn and thus restore the status of tin(p) as 
the distribution function. 



APPENDIX B. MISCELLANEOUS FORMULAE. 

1. If we adopt the Boltzmann shape of the distribution function, then the product of the statistical weights 
n(pi)n(p 2 ) becomes a function of p — p\ + p 2 and it is possible to integrate / = p\ — p 2 out. We have an integral, 

J d i Pl d i p 2 f(p)S + (p 2 1 - m 2 )5 + (pj -m 2 )= l -J d^pfip) J dHS(p 2 + I 2 - Am 2 )5{pl) 




(B.l) 



where in the second equation we integrate using a special reference frame where p = 0, and restore the invariant form 
in the final answer. 

2. The correlation functions of Sec. are expressed via the emission function I{k,N,XNM) which is computed 
below: 
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I(k,N,X NM ) = / d 4 qS+(q 2 - m 2 )e-^ e iqx Jl 



4m 2 



(B.2) 



where (3 — /3jy, u — un, and x — xnm- The integral is convenient to compute in a special reference frame 1Z where 

o o 

the time-like vector k has no spatial components, k = 0, ko — m. Using the delta-function to integrate out q , we 
arrive at 



2;r 



I(k,N,XNAi) = / dcf) I sinada 



qdq 



o 2^/q 2 + m A 

x exp{— \J q 2 + m 2 (P — ixo) + <z[/?|w| cosa — i \ x \ (cos a cos ip — sin a sin tp cos (j))]} 



(B.3) 



The first step is to integrate over <fi which results in the Bessel function, 27rJo(|£c|<7sinasin , (/>). The next integration 
follows a known formula, 



da sin a e 



ib cos a 



Jo(asma) = 2 



sin Va 2 + b 2 
Va 2 + b 2 



(B.4) 



which holds for arbitrary complex a and b. In our case, a 2 + b 2 = —{(3u — ix) 2 = —U 2 where we have introduced the 
complex four- vector — f3u^ — ix^ . Finally, changing the variable of integration, q = msinhw, we arrive at 



I(k,N,x NM ) = 



2nm 



2nm 



\U\ Jo 

dv (coshxcoshw- i) e - m ^ COBh » 



(coshv - l)dv e" mC/nCOshl, sinh(m|Lr| sinhw) 

dv (coshw - l) e -"»v^coBh(t;+x) 



v U 2 sinh x 

o 

The Lorentz-invariant expressions for the quantities defined in the reference frame 1Z are 



iUq = (kU) = mvU 2 coshx, m\U\ = \J (kll) 2 — m 2 U 2 = mvU 2 sinh \ 



(B.5) 



(B.6) 



The first of the integrals in Eq. (B.5) is calculated exactly. The second one can be estimated by means of integration 
by parts. Since the integrand has the root of the second order at v = 0, we have to integrate by parts three times 
before the first non-vanishing term shows up. The net yield of this procedure is 

47T77J 2 



Airm 



\U\ Jo 



dv Sinh 2 ^-'"V^cosh^+x) 

2 



-mUn 



m 4 \U\ 4 



(B.7) 



Since mC/o = (kU) = (ku)/T — i(fcxjvM); this term is strongly suppressed (both by the exponent and by at least the 
fourth power of the small number m/kt). Finally, we have to high accuracy, 



I(k,N, x NM ) 



2nm 



coth xKi (mVu 2 ) 



1 



sinh x 



KoimVu 2 ) 



(B. 



Here, K n (z) is the standard notation for the modified Bessel function, and we remind the reader that U 2 and x 

o o 

(as well as Uq and \U\) all are complex quantities. When x = xnm — 0, we have U 2 — (3 , coshx — (ku)/m — 
(rnt/m) cosh(6> — ry) 3> 1, and consequently, sinhx ^ 1, cothx ~ 1- Therefore, 



I(k,N,0) = 2irm 



2 Ki (J3m) 



(3m 



(B.9) 



In the particular cases considered in Sec. 0, the second term in Eq. (B^S) is always smaller than the first one by the 
factor m/k t . An important property of the final answer (B.8) is that the quantity 
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U 2 = (3 2 ( Xl ) - ( Xl - x 2 f - 2i0(x 1 )(u(x 1 ) ■ ( Xl - x 2 )) 



in the argument of function K\, does not depend on the large parameter k t /m. Therefore, the saddle point of the 
integration over the coordinate rapidities r\ is never affected by this function. 
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